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It has been noted for a long time that, in some circumstances, massless 
neutrinos may be mixed in the leptonic charged current. Conventional 
1 J ■ neutrino oscillation searches in vacuum are insensitive to this mixing. We 
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H ' discuss the effects of resonant massless-neutrino conversions in the dense 

medium of a supernova. In particular, we show how the detected z/g energy 
spectra from SN1987a and the supernova r-process nucleosynthesis may 
be used to provide very stringent constraints on the mixing of massless 
neutrinos. 
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1 Introduction 



In the original scenario developed by Mikheyev and Smirnov the resonant neu- 
trino conversion in matter requires non-degenerate neutrino masses and non-vanishing 
mixing angles in vacuum. In the basis of two neutrino flavour eigenstates, the evolution 
Hamiltonian describing the neutrino propagation in matter is given by 

^ = ( 5' ) ' " = (1) 

Hea = Hae = 77^ slu 29, 

where Vg and are the well-known Wolfenstein diagonal matter potentials arising from 
coherent neutrino scatterings off matter particles p[. In Eq. (1), E,y is the neutrino 
energy, and and 6 are the neutrino mass-squared difference and the mixing angle 
in vacuum, respectively. One can see that the effective mixing in matter between z/g 
and Ua states is induced by the "vacuum" term (drn^/AE^) sin 26*. 

It has been noticed for a long time that the presence of SU{2) ^ U{1) isosinglet 
neutral heavy leptons in general leads to flavour-changing neutral-current (FCNC) 
interactions of neutrinos [^. As a result, there can be non-trivial leptonic mixing (and 
CP violation) involving the conventional isodoublet neutrinos even in models where 
these neutrinos remain strictly massless, as in the Standard Model, due to an exactly 
conserved lepton number . The non- vanishing massless-neutrino mixing angles arise 
due to the presence of extra heavy gauge singlet neutral states. In this scenario the 
interaction of massless neutrinos with matter constituents gives rise to a non-trivial 
neutrino evolution Hamiltonian, analogous to Eq. (|l|) P, ^, which can mix the neutrino 
identities p, ^. This Hamiltonian is characterized by a new type of weak potentials 
whose diagonal and off-diagonal matrix elements will be discussed below. 

The implications of both standard and non-standard neutrino interactions for 



the neutrino propagation in dense media have been extensively studied |]I0|, [111. In 
particular, the birth of neutrino astronomy, with the detection of neutrinos from the 
Sun ||12[ and SN1987a has offered the opportunity to probe various neutrino 
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properties, such as neutrino masses and mixings, neutrino lifetimes, neutrino magnetic 
moments, and generically, any non-standard interactions of neutrinos. 

In this paper we focus on the particular scenario of massless-neutrino mixing 
suggested in Ref. [||. This scenario can be relevant only for the neutrino propaga- 
tion in strongly-neutronized media. Such media exist perhaps only in supernovae. 
We show how to probe the mixing in the light neutrino sector by considering two 
different aspects of the supernova process. We examine how the massless-neutrino 
conversion of the type z/g ^ can affect the detected z/g energy spectra [0, |I6|, |1^ 
from SN1987a. We also consider the implications of such conversions for the super- 
nova r-process nucleosynthesis, following the same lines of reasoning adopted in Ref. 
||18|| . Rather stringent limits on the mixing of massless neutrinos may be derived from 
both considerations. These limits are very remarkable because the mixing of massless 
neutrinos cannot be sharply constrained through neutrino oscillation searches. Being 
strictly massless, these neutrinos cannot develop any phase difference in vacuum, and 
as a result, neutrino oscillations can not occur. 

In Sec. 2 we give a quick reminder on the theoretical framework of Ref. 0. In 
Sec. 3 we present the general features of the resonant massless-neutrino conversion in 
matter. Sec. 4 discusses resonant conversions of massless-neutrinos in supernovae and 
the implications of such conversions for supernova neutrino detection and r-process 
nucleosynthesis. We summarize our results and conclude in Sec. 5. 



2 The Massless and Mixed Neutrino Model 



In the Standard Model the absence of right-handed neutrino states naturally 
implies that neutrinos stay massless to all orders of perturbation even after the gauge 
symmetry breaking and there are no Cabibbo-Kobayashi-Maskawa-like |]T^ mixing ma- 
trices in the weak leptonic charged current. In this case the total lepton number L 



comes out as an accidental symmetry |2^ due to the gauge structure and renormaliz- 
ability of the theory. 

On the other hand, any number of gauge singlet neutral leptons can be introduced 
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since they do not carry triangle anomaly These extra states can arise in left-right 
symmetric, grand-unified or superstring-inspired models @, ^ |2l], In this case 

the lepton number is no longer an accidental symmetry and it may be imposed by 
hand. The simplest such scheme [§], §|, 01 contains three two-component gauge singlet 
neutral leptons 5* added to the three right-handed neutrino components u'^ present in 
SO(IO). For definiteness we consider this model at the SU{2)(^U{1) level. The assumed 
conservation of lepton number leads to a neutral mass matrix with the following texture 
in the basis (z/, z/'^, S): 

D \ 

D'^ M , (2) 

; 

where the Dirac matrix D describes the coupling between the weak doublet v and 
the singlet v'^, and where the other Dirac matrix M connects the singlet states v'^ 
and S. It is easy to see that, as expected, the three conventional neutrinos remain 
massless, while the other six neutral 2-component leptons combine into three heavy 
Dirac fermions 0, |^. 

The phenomenological implications of this picture are manifest when considering 
the resulting charged-current (CC) Lagrangian in the massless-neutrino sector: 

Ccc = -^W^eaLl^,Kail'iL + ^■c., a = e,fi,T, i = 1,2,3, (3) 

where the mixing matrix K is not unitary, since it is a sub-matrix of the full rectangular 
matrix including also the heavy states [Q. Therefore, the non-diagonal elements of 
the matrix K cannot be rotated away through a redefinition of the massless-neutrino 
fields. In this way a non-vanishing mixing arises among the massless neutrinos. The 
corresponding form of the neutral-current (NC) Lagrangian for the massless-neutrino 
sector is 

= o^^f/j Zf,PijUiLjf,PjL , (4) 

z cos 

where P = K'^K. Unlike in the Standard Model, the matrix P is diagonal but 
generation-dependent, signalling the violation of weak universality. 

For definiteness, we later on use an explicit parametrization of the matrix K, 
confining ourselves to the case of two (massless) neutrinos u. We may write the mixing 
matrix i^' as @, ||] 

K = RAf, (5) 
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where i? is a 2 x 2 rotation matrix. 



/ cos 6 sin 6* \ 

^= • . J' 
\ — sm 6 cos U ) 

and where the diagonal matrix, 

describes the effective non-orthogonahty of the two neutrino flavours, i.e., 

(t'elt'/^.r) = — sin cos0(A/'f — ■^2,z)- ^hc corresponding NC couplings in Eq. (H) are 

now expressed through 

P=M^. (8) 

It is also convenient to define 

M^^{l + h^)-\ z = 1,2(3), (9) 

where the hi parameters reflect the deviation from the standard neutrino coupling. 

Before entering into the discussion of the resonant massless-neutrino conversion 
in matter, we describe the present upper limits from laboratory experiments on the 
relevant parameters hf and 9. We first note that the laboratory limits on the leptonic 
mixing angle 9 are rather weak since no oscillations between two strictly massless 
neutrinos can develop in vacuum. However, although not strictly justified a priori 
from the point of view of laboratory constraints, we will assume the small-mixing 
angle approximation. This will be justified a posteriori in view of our results. In this 
way we have Ui ~ z/q [a = e,/i(r)], and we can analogously interpret hf as h^. 

There have been extensive studies of experimental universality tests which restrict 
the parameters h^. For the case of h"^ one can still allow values in the range of a few 
percent |P, whereas the constraints on and /i^ are more stringent. Therefore, 



from now on we focus on the {ue, Vt) system, for which the universality limits are 
the weakest. Moreover, the present experimental situation cannot exclude that the 
difference h'^ — hi can be positive as required later on in our discussion. 



3 Resonant Massless-Neutrino Conversion 
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Here we briefly recall the main features of the resonant conversion mechanism of 
massless neutrinos emerging from the previous scenario. For convenience we choose to 
write the system of Schroedinger equations, which describe the propagation of the two 
neutrinos in matter, in the basis defined as ^ 

z>, = [RM-^R'U Vh. a,b = e, r. (10) 

Although this basis is somehow artificial, it almost coincides with the flavour basis 
for small lepton universality violation or small mixing angle 6. In this basis, the 
Schroedinger equations can be written as 

dr \ aJ ^mw \h„ hJ\ AJ' 

where ,- are the amplitudes corresponding to the neutrino states in the basis of Eq. 
(p!0|), Gj? is the Fermi constant, p is the matter density, and m^r is the nucleon mass. 
The entries of the evolution Hamiltonian are now given by 

He = YeiKc'+Ks'f-lYr^iAfy+Afls'), 

Her = Hre=[YeiKc''+Ks^)-^Yr,iK+K)]iK-K)sC, (12) 

Hr = FeA^(AC-AC^)-^F„(ACV+ACV), 

where for brevity we have used the short-hand notations s = sin 6' and c = cos^. 
In an electrically neutral medium. Ye and Yn are defined as 

rj 

Ye = ^ , y„ = 1 - , (13) 

Ue + Un 

where ne and the net electron and the neutron number densities in matter, 

respectively. Note that the evolution matrix has no energy dependence, which implies 
that for the corresponding antineutrino system (i/ej^r) this matrix just changes its 
overall sign. Clearly, in this scenario, resonant neutrino conversion can also occur 
provided the condition He = Hr is fulfilled In fact, the same resonance condition 
holds for both z/g ^ and Ue ^ i^r channels. As a result, in the thermal phase 
of supernova neutrino emission, both neutrinos and antineutrinos can simultaneously 
undergo this resonance. This will be very important for our subsequent discussion in 
Sec. 4. 



5 



In order to simplify Eq. (0), we take advantage of the small parameters h"^ 
expected from the universality constraints. With the previous assumption of small 6, 
we obtain 



I— I ~ I = \/2Gf I 1 1 o 

dr\Aj mM\\r,{Y^-Y,)sin2e -\Y^{1 - hi) 



(14) 



where the parameter rj is defined as 



1 



(15) 



The mixing angle 6'^ and the neutrino oscillation length in matter are given by 

7]\Yn~Y,fsin^2e 



''''' ^^"^ {Ye - vYnY + v'iYn ' Y^f siu^ 29' 

27r 

^Gpip/mNW, - r^Y,^Y + r/2(F„ - Y^f sin^ 29]^/^ ' 



(16) 
(17) 



respectively. 

The resonance condition now reads 



Ye = VYn. 



Here we should stress that a positive value of t] is necessary for the above equation to 
hold. Moreover, due to the bounds on the lepton universality violation, rj < O(10~^), 
the condition in Eq. ([18D can be fulfilled only in a strongly-neutronized medium. 
This is why the present mechanism cannot work in the matter background of the Sun 
{Yn < 0.33) 1^, or Earth (y„ ~ 0.5). On the other hand, the material composition 
just above the neutrinosphere in type II supernovae (Fg ^ F„) can satisfy Eq. (0), as 
shown later. 

In our subsequent discussion, we will employ the simple Landau-Zener approx- 
imation |25| to estimate the conversion probability after the neutrinos cross the 
resonance. Under this approximation, the probability for <-> Uj. and ^ Vj- conver- 
sions is given by 



1 



exp 



2 Ll 



exp 



-32 X 



5r 



2Hsm29, H 



10i2g/cmV 
dlnK 



10- 



—\ smH9 
cmj 



dr 



(19) 
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where is the neutrino oscillation length at resonance and p^es are the corresponding 
matter density. In deriving the above equation, we have used y„ ~ 1 for the neutron 
abundance near resonance. Notice that for 5r/U^ > 1 resonant neutrino conversion 
will be adiabatic It is also important to note that the conversion probability does 
not depend on the neutrino energy [cf. Eqs. (|TT|) and (|T^]. 

4 Massless-Neutrino Conversion in Supernovae 
4.1 Neutrino Emission and Y^, Profile in Supernovae 

A supernova occurs when the core of a massive star collapses into a compact 
neutron star. Almost all of the gravitational binding energy of the final neutron star 
is radiated in Ue, i^e, J^fi, i^fi, J^r, and Ur- The last four neutrino species are created by 
thermal pair production processes inside the neutron star. On the other hand, although 
most of the and z/g are produced in pairs, there is a net excess of z/g over u^, due to 
the neutronization and deleptonization of the core through + p ^ n + Ue. Because 
all these neutrinos have intense neutral-current scatterings on the free nucleons inside 
the neutron star, the net lepton number carried by z/g can escape from the neutron star 
only through diffusion. Therefore, we expect to see the strongest deleptonization effect 
near the neutrinosphere, where neutrinos stop diffusing and begin free-streaming. 

We can estimate the electron fraction near the neutrinosphere as follows. From 
the approximate chemical equilibrium for e~, p, n, and z/g, we have 

/ie- + yUp ~ /in, (20) 

where for example, /ig- is the electron chemical potential, and where we have set 
/ij,^ ~ 0. For non-relativistic nucleons, we can write 

- - exp ^—^j , (21) 

where T is the temperature, and where we have neglected the neutron-proton mass 

difference. The electron fraction is then given by 

n„ 1 . . 



The chemical potential for relativistic and degenerate electrons near the neutrinosphere 
is approximately given by 

/ie- ~ (?>Ti^n,Y'^ ^ 51.6(FePi2)'/' MeV, (23) 

where pi2 is the matter density in units of 10^^ g cm~^. For typical conditions near 
the neutrinosphere, T ~ 4 MeV and pi2 ~ 10, by solving Eqs. (^) and (^), we find 
~ 6 X 10~^, in good agreement with the numerical supernova models. Therefore, we 
can expect resonant massless-neutrino conversions to occur above the neutrinosphere 
as long as the lepton non-universality parameter > 6 x 10~^ [cf. Eq. (|T^)]. 

Above the neutrinosphere, the approximate chemical equilibrium between i/g and 
matter no longer holds. The electron fraction is determined by the following reactions: 

Pe + n ^ P + e", (24) 
Ve+V ^ n + e"*". (25) 



In fact, Qian et al. [18| have shown that Yg above the neutrinosphere is given by 



where for example, Aj,^^ is the rate for the forward reaction in Eq. (p^. In particular, 
because X^-p and Ae+„ quickly decrease with the temperature, the asymptotic value of 
Ye at large radii is approximately given by 

Ye ^ , . (27) 

Therefore, the asymptotic electron fraction above the neutrinosphere is essentially de- 
termined by the characteristics of the z/g and z/g fluxes, such as their luminosities and 
energy distributions. 

The individual neutrino luminosities in supernovae are approximately the same: 

However, the individual neutrino energy distributions are very different. This is be- 
cause these neutrinos have different abilities to exchange energy with the neutron star 
material, and thermally decouple at different temperatures inside the neutron star. 
Unlike z/g and z/g, z/,-(^) and z/^(^) are not energetic enough to have charged-current ab- 
sorptions on the free nucleons inside the neutron star. Furthermore, between z/g and z/g. 
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Vf. have more frequent absorptions due to the high neutron abundance in the neutron 
star matter. As a result, z/t-(^) and z/t-(^) thermally decouple at the highest temperature, 
and Ve decouple at the lowest temperature. Correspondingly, the average neutrino 
energies satisfy the following hierarchy: 

< (Ep,) < l,E,^J ^ {E,^J. (29) 

Typically, the average supernova neutrino energies are: 

^ 11 MeV, (EpJ ^ 16 MeV, {E^^J ^ {E,^J ^ 25 MeV. (30) 

Now we can understand the electron fraction profile in supernovae as illustrated 
in Fig. 1. In this figure, we plot the typical electron fraction and density profiles in 
Wilson's supernova model at time t > 1 s after the explosion. The solid line is for the 
density, and the dotted line is for Yf.. As we can see, the minimum value of occurs 
near the neutrinosphere. Above the neutrinosphere, the electron fraction is set by the 
reactions in Eqs. (0) and (p5|). At large radii, it reaches an asymptotic value much 
larger than the minimum Yg. 



From the above discussion of neutrino emission and Y^ profile in supernovae, we 
find that it is interesting to study massless-neutrino conversion in supernovae. First 
of all, the resonance condition for such conversion, Eq. (ITH) can be fulfilled above the 
neutrinosphere for t] ~ 0.01. Furthermore, conversion between (z/,-) and z/g (z/g) can 
alter the supernova neutrino characteristics, especially the average neutrino energies 
in Eq. (pOl). We can gauge the potential to use supernovae as a sensitive probe of the 
mixing between massless neutrinos by estimating the adiabatic condition for resonant 
massless-neutrino conversion. For rj ~ 10~^, the resonances occur at densities p ~ 10^^- 
]^Qi3g Q]2i-3^ just above the neutrinosphere. The corresponding scale height for Y^ is 
H ~ 1-10 km. From Eq. ([T9|), we see that massless neutrinos can be adiabatically 
converted for sin^ 29 > 10^^-10"^. In the following subsections, we discuss two possible 
ways to probe the mixing between massless neutrinos in supernovae. 



4.2 Detection of z>e from SN1987a 
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The Kamiokande II and 1MB detectors observed 11 and 8 Ve events, respectively, 
from SN1987a [9,10]. An estimate of the average supernova energy can be made 
from the detection data, although the obtained estimate should be taken with caution, 
considering the poor statistics and the marginal agreement between the two sets of 
data. Nevertheless, if we adopt the standard average neutrino energies predicted by 
the numerical supernova models, then a significant amount of conversion between Vj. 
and Vf. can probably be ruled out. This is because the average Ve energy inferred 
from the detection data is much smaller than the average z/^ energy predicted by the 
numerical supernova models. Specifically, with conversion, the z/g flux at the 

detectors would be given by 

= 0°^(1-P) + 0°^P, (31) 

where (pp^ and 0°^ are the z/g and Ur fluxes in the absence of neutrino conversions, 
respectively, and P is the conversion probability. For large P, based on predictions 
from numerical supernova models, the z/g energy spectra at the detectors would have 
been significantly harder than detected in the case of SN1987a. From the detection 
data, Smirnov et al. argued that the probability for ^ ^t(p) conversion should 
be less than 0.35. 

We can apply the same argument to constrain the mixing between massless neu- 
trinos. Using the density and profiles from Wilson's supernova model in Fig. 1, we 
plot in Fig. 2 two contours of the conversion probability in the (77, sin^ 26) parameter 
space. The solid line is for a conversion probability of P ~ 0.5, and the dashed line 
is for P 0.35. We can conclude that mixing between massless z/g (z/g) and v^- {vr) 
at a level of sin 29 > IQ-^ is ruled out for t] > IQ-^ due to the non-observation of 
unexpectedly hard energy spectra from SN1987a. Such a stringent upper limit on 
the mixing angle 6 justifies the approximation we have made in deriving Eq. (0). 



4.3 Supernova r-process Nucleosynthesis 



Now we consider the effect of massless-neutrino conversions on the the supernova 
r-process nucleosynthesis. The r-process is responsible for synthesizing about a half 
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of the heavy elements with mass number A > 70 in nature. It has been proposed 
that the r-process occurs in the region above the neutrinosphere in supernovae when 
significant neutrino fluxes are still coming from the neutron star ^6]. A necessary 



condition required for the r-process is Y^. < 0.5 in the nucleosynthesis region. As we 
have discussed previously, the Y^ value at large radii above the neutrinosphere, where 
the r-process nucleosynthesis takes place, is determined by the neutrino absorption 
rates Aj,^„ and Ap^p. In turn, these rates depend on the Ue and luminosities and 
energy distributions. 

Qualitatively, we can argue that these rates are proportional to the product of 
the neutrino luminosity and average neutrino energy. This is because the neutrino 
absorption rate is given by 

XuN - <PA^un) ^ ^{ED ^ L,{E,) , (32) 

where (p^, is the neutrino flux, (7^]\r oc is the neutrino absorption cross section, and 
( ) denotes the averaging over the neutrino energy distribution. Therefore, the Y^ in 
the nucleosynthesis region is approximately given by 

Using the average energies in Eq. (0), we obtain Y^. ^ 0.41, in good agreement with 
the numerical supernova models. 

However, in the presence of massless-neutrino conversion, average energies of both 
Ue and Ug can be affected. The corresponding Yg in the nucleosynthesis region is given 
by ^ 

" l + {E,X,/{EuX,' ^^^^ 



where 



{E,J^^^{E,J{1-P) + {E,^)P, (35) 
{EA^^^{EAil-P) + {EAP. 

Due to the the simultaneous occurrence of resonant z/g ^ Ur and ^ i^r conversions, 
there is a trend to equalize the average z/g and Ue energies, and as a result, to increase Ye 
with respect to the case with no neutrino or antineutrino conversions. For conversion 
probabilities of P ^ 0.15, 0.3, and 0.8, we obtain Ye ^ 0.43, 0.45, and 0.49. In Fig. 3, 
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we present the contour lines corresponding to these Ye values using the density and 
profiles in Wilson's supernova model. The dotted, dashed, and solid lines in this figure 
are for Y^ ^ 0.43, 0.45, and 0.49, respectively. 

In order for any r-process nucleosynthesis to occur, the Y^ in the nucleosynthesis 
region must be less than 0.5. However, in the most recent r-process model by Woosley 



et al. [^, many of the r-process nuclei are produced only for Y^ < 0.45. If we take 
Ye < 0.45 as a criterion for a successful r-process, then mixing between Vf. {ug) and 
z/^ (z7^) at a level of sin^ 29 > 10^^ is excluded for t] > 10^^. This excluded region is 
similar to the previous one from considering the detection of Ue from SN1987a, because 
the limits on the conversion probability are about the same in both cases. However, 
we note that if the r-process indeed occurs in supernovae, then the consequent limits 
on the mixing between massless neutrinos are much less dependent on the predicted 
average neutrino energies than the previous limits obtained by considering the Ue energy 
spectra from SN1987a. This is because the r-process argument relies only on the ratio 
of the average neutrino energies [cf. Eq. (|33D ]. 

4.4 Comparison of MSW and Massless-Neutrino Conversion 
Mechanisms in Supernovae 



It is instructive at this stage to compare the effects of resonant massless-neutrino 
conversions with those of the standard MSW mechanism in supernovae. To simplify 
this comparison, we will assume small vacuum mixing angles (^^ ^ 1) in both cases. 
We first note that in the MSW scenario IQ, for a given sign of 6m^ (e.g., > for 
■m,^^ > m^J, only one kind of resonant conversion, either z/g ^ Ur (for ^^"^ > 0), or 
z7g (for < 0), can occur. If Sm^ > the MSW mechanism would not alter 

the z/g energy spectra from SN1987a, and therefore, no constraints on neutrino masses 
and mixings can be obtained for this mechanism from the detection data (assuming all 
the events were due to z/g). In contrast, severe constraints on massive-neutrino mixing 
can be obtained in this case by requiring Ye < 0.5 in the nucleosynthesis region to 
allow a successful r-process On the other hand, if < the MSW mechanism 
could significantly modify the z/g energy spectra and generate an excess of energetic 
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z/g. As a result, the parameter region which would give large probabilities for v^. Ur 
conversion can be possibly excluded by combining the predicted average supernova 
neutrino energies and the SN1987a detection data. On the contrary, significant Ue ^ i^r 
conversion would tend to decrease Y^, in the nucleosynthesis region (see eq. (0)) and 
therefore, would not confiict with the supernova r-process nucleosynthesis scenario [25]. 
As we can see, one can only use either the SN1987a detection data (for < 0), or 
the supernova r-process nucleosynthesis (for Sm'^ > 0) to constrain neutrino masses 
and mixings in the MSW mechanism. 

In contrast, in the case of massless-neutrino conversions, we have seen that for 
rj > 10^^, both Ue ^ v-r and z/g ^ z^r conversions can occur in supernovae. Therefore, 
both the SN1987a detection data and the supernova r-process nucleosynthesis should 
be considered in order to constrain the mixing of massless neutrinos. Of course, if 
< or <^ 10^^, then no resonant massless neutrino conversions would occur in 
supernovae. The constraints on massless-neutrino mixing in this case are perhaps hard 
to obtain by any means. 

It is interesting to note that simultaneous z/g ^ and z/g Vt conversions 
would give rise to distinctive supernova neutrino signals in large volume detectors, 
such as super-Kamiokande [26], SNO [27], and LVD [28]. For example, in the super- 
Kamiokande detector, the energy distributions for both the isotropic z>e events and 
the forward-peaked z/g events would be altered. With enough statistics, such detectors 
may be able to distinguish the massless-neutrino conversion scenario from the standard 
MSW mechanism, should any neutrino conversion indeed occur in supernovae. 



5 Conclusions 



We have discussed the possibility of probing the mixing between massless neutrinos 
described in the theoretical scheme in Sec. 2. Due to the relatively stringent labo- 
ratory bounds on the weak universality violation, the supernova matter background 
seems to be the unique site where resonant conversions of massless neutrinos can take 
place. By considering the detection of z>g from SN1987a and the supernova r-process 
nucleosynthesis, we have obtained stringent limits on the mixing between massless Ve 
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(z/g) and v-r {pr) presented in Figs. 2 and 3. These limits, at a level of sin^2^^ < 10~^, 
are rather remarkable, because the usual laboratory methods to constrain neutrino 
mixing through vacuum neutrino oscillation searches are totally insensitive to the mix- 
ing between massless neutrinos. Indeed, the supernova limits we have obtained for the 
mixing between massless (Pg) and (P,-) are orders of magnitude more stringent 
than the typical limits on massive-neutrino mixing from laboratory neutrino oscillation 
searches. 

Finally, we hope that our discussions of massless-neutrino conversions in super- 
novae serve to highlight the interest in sharpening the laboratory limits on universality 
violation and/or pinning down more accurate supernova models. 
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Fig. 1 
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Figure 1: Typical matter density (solid line) and (dotted line) profiles in Wilson's 
numerical supernova model at i > 1 s after the explosion. 
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Fig. 2 
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Figure 2: Constraints on massless-neutrino mixing from the detected SN1987a Ve 
energy spectra. The region to the right of the dashed (sohd) hnes are excluded by the 
detection data for an allowed conversion probability of P < 0.35 (0.5). 
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Fig. 3 
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Figure 3: Constraints on massless-neutrino mixing from the supernova r-process nu- 
cleosynthesis. The region to the right of the dotted, dashed and sohd hues are exclued 
for the required values of Ye < 0.43, 0.45, and 0.49, respectively, in the r-process. 
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